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ABSTRACT 

RS Ophiuchi is a recurrent nova system that experiences outbursts every ^20 years, implying 
accretion at a high rate onto a massive white dwarf. However, previous X-ray observations of the 
system in quiescence have detected only faint emission that is difficult to reconcile with the high 
accretion rate predicted by nova theory for such frequent outbursts. Here, we use new Chandra and 
XMM-Newton observations obtained 537 and 744 days after the 2006 outburst to constrain both the 
accretion rate onto the white dwarf and the properties of the nova ejecta at these times. We detect 
low level UV variability with the XMM-Newton Optical Monitor on day 744 that is consistent with 
accretion disk flickering, and use this to place a lower limit on the accretion rate. The X-ray spectra in 
both observations are well described by a two component thermal plasma model. The first component 
originates in the nova shell, which can emit X-rays for up to a decade after the outburst. The other 
component likely arises in the accretion disk boundary layer, and can be equally well fit by a single 
temperature plasma or a cooling flow model. Although the flux of the single temperature model implies 
an accretion rate that is 40 times lower than theoretical predictions for RS Oph, the best fit cooling 
flow model implies M < 1.2 x 10~ 8 M yr _1 537 days after the outburst, which is within a factor of 2 
of the theoretical accretion rate required to power an outburst every 20 years. Furthermore, we place 
an upper limit on the accretion rate through an optically thick region of the boundary layer of 2.0 x 
10~ 8 M yr . Thus, these new quiescence data are consistent with the accretion rate expectations 
of nova theory. Finally, we discuss the possible origins of the low temperature associated with the 
accretion component, which is a factor of 10 lower than in T CrB, an otherwise similar recurrent nova. 
Subject headings: stars: white dwarfs, X-rays: stars 



> 
OS 

SO 
in 

o 



1. INTRODUCTION 

Recurrent novae are a rare subclass of cataclysmic vari- 
able stars (CVs; white dwarfs accreting material from 
a binary companion) in which more than one classical 
nova-type outburst has been observed. These outbursts 
are believed to be due to a thermonuclear runaway on 
the surface of the white dwarf, which commences once 
a critical pressure is reached at the base of the shell 
of accreted material. All classical novae are expected 
to recur on timescales from 100,000 years to just a few 
decades. The most important physical parameters con- 
trolling this recurrence timescale are the white dwarf 
mass, an d the mass accret ion rate from the secondary 
(see e.g. lYaron et alJl200l . RS Oph is a particularly 
prolific recurrent nova, with detected outbursts in 1898, 
1907 , 1933, 1945, 1958, 1967, 1985 and 2006 (ISchaeferl 
2010( 1. The short time between outbursts (~20 yrs) sug- 
gests that RS Oph hosts a massive white dwarf accreting 
material at a high rate. 



thomas.nelson@nasa.gov 

1 ORESST and X-ray Astrophysics Laboratory NASA/GSFC, 
Greenbelt, MD 20771, USA 

2 Center for Space Science and Technology, University of 
Maryland Baltimore County, 1000 Hilltop Circle, Baltimore MD 
21250, USA 

3 INAF - Osservatorio Astronomico di Padova, vicolo Osser- 
vatorio, 5, 1-35122 Padova, Italy 

4 Department of Astronomy, University of Wisconsin- 
Madison, 475 N Charter St, Madison, WI, 53706, USA 

5 Harvard-Smithsonian Center for Astrophysics, 60 Garden St. 
MS15, Cambridge, MS 02138, USA 

6 Columbia Astrophysics Laboratory, 550 W. 220th St, 1027 
Pupin Hall, Columbia University, New York, NY 10027, USA 



Observations of RS Oph during its 2006 outburst al- 
lowed direct observational constraints to be pla ced on the 
mass of the white dwarf. iNelson et~a l. (2008) presented 
high resolution X-ray spectra obtained during the super- 
soft phase (the phase where the still burning hydrogen 
shell is directly observable) , and found a very high effec- 
tive temperature for the source (~800,000 K), consistent 
with shell burning on a white dwarf of at leas t 1.2 M 
(see e.g. lPaczvnski|[l97lt iNomoto et all 120071 ). Subse- 
quent an alysis of the Swift data obtained during the same 
epoch bv lOsborne et alj (|2011f ) led those authors to also 
conclude that the whit e dwarf mass must be more than 
1.2 M Q . Furthermore. iSokoloski et al.l (|2006[ ) estimated 
the mass of the ejected envelope based on the hard X-ray 
light curve of the outburst, and found that no more than 
a few 10 -7 M could have been ejected. This is con- 
sistent with the predicted ejecta ma ss for a nova even t 
on a 1.4 M WD (-2 x 10~ 7 M , lYaron et all 120051) . 
Thus, the white dwarf in RS Oph is extremely massive, 
in agreement with the first expectation for systems with 
frequent nova outbursts. 

However, there is less evidence of the high accretion 
rate required to power an outburst every ~20 yrs. Cat- 
aclysmic variables are well known X-ray sources, and 
are understood to be powered by accretion. For a 
slowly rotating white dwarf, approximately half of the 
available accretion luminosity is emitted in the accre- 
tion disk boundary layer at X-ray wavelengths (see e. 
IShakura fc Sunvaevl[l97l ILvnden-Bell fc Prinrid[l974 
Therefore, a measurement of the X-ray luminosity can 
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be used to estimate the accretion rate as follows: 



where M is the mass accretion rate, Lx is the X-ray lu- 
minosity, and Rwd and Mwd and the radius and mass of 
the white dwarf, respectively. The same method c an be 
applie d to existing X-ray obs ervations of RS Oph. 10 rid 
(fl99l and IQrio et all (f200l presented ROSAT PSPC 
observations of RS Oph carried out 6 and 7 years after 
the 1985 outburst, respectively. In those observations, 
the system was detected as rather faint, soft X-ray source 
(all photons with energies <1.5 keV), which varied in 
count rate between the two observations by a factor of 
3. The derived luminosities in the ROSAT energy range 
(0.2-2.4 keV), based on best fitting thermal plasma mod- 
els and assuming a distance to the system of 1.6 kpc, 
were 3 x 10 31 to 5 x 10 32 erg s -1 . Using Equation 1, 
and assumi ng a white dwarf ma ss of 1.3 Mq and radius 3 
x 10 8 cm (jAlthaus et al.l 120051 ) . this low luminosity im- 
plies a mass accretion rate in the range 10 -12 to a few 
KT 11 M Q yr" 1 . A similar range of luminositie s, and 
hence accretion rates, was found bv lMukail (2009) using 
a 35 ks ASCA observation of the system carried out in 
1997 (~12 yrs after the 1985 outburst). Such low rates 
are completely at odds with the short nova outburst re- 
currence time (~20 years), which requires M >10~ 8 M Q 
yr" 1 for a 1.4 M Q WD, and >10~ 7 M w yr" 1 if the mass 
of the WD is only 1.25 M (jYaron et al.|[2005T) . 

RS Oph is also categorized as a symbiotic star i.e. a 
star whose binary nature is revealed through the pres- 
ence of two distinct components in its spectrum, one hot 
and one cool. Most of these systems are comprised of a 
white dwarf accreting material from a red giant compan- 
ion in a wide binary, with orbital periods in the range of 
a few years to a few decades. Symbiotic stars can there- 
fore be considered larger scale cousins of the cataclysmic 
variables. The mode of accretion onto the white dwarf 
has not been definitively identified in symbiotic stars. 
A small number of systems, including RS Oph, exhibit 
strong optical flickering (stochastic brightness fluctua- 
tions on timescales of mi nutes to hours) , whic h is also 
seen in CVs with disks (jSokoloski et al.ll200l . These 
brightness fluctuations are thought to be due to the tur- 
bulent motio n of material as it passes through the disk 
(|Bruchlll99l . and their presence is taken as evidenc e of 
an accretion disk in those systems. [Sokoloski (200l|) re- 
ported the presence of minute time scale variability in 
the He II A4686 line in RS Oph in quiescence. The high 
ionization potential of this ion (56 eV) requires a source 
of soft X-rays, most likely the accretion disk boundary 
layer. Addition ally, the launching of r adio jets during the 
2006 outburst (jSokoloski et alj|2008t iRup'en et al.ll2008D 
can also be considered evidence of the presence of an ac- 
cretion disk, since almost all mechanisms for producing 
jets require a disk to anchor the magnetic fields that col - 
limate the outflow. We refer the reader to IWvnnl (|2008D 
for a detailed discussion of the possible modes of accre- 
tion in RS Oph. 

Although the presence of a disk in RS Oph seems 
likely, the mechanism by which material is fed into it 
is unknown. Most of the giants in symbiotic stars do 
not appear to fill their Roche Lobes, although a number 



of systems show near-infrared ellipsoidal variations, sug- 
gesting that they ar e at least partially tidally distorted 
(|Mikolaiewskal I2007D . If Roche lobe overflow (RLO) is 
not occurring, it is normally assumed that material is 
fed to the white dwarf through Bondi-Hoyle accretion of 
the red giant wind. Since this mode of accretion is much 
less efficient than RLO, it is unclear that it can provide 
enough mass to produce an outburst in RS Oph every 20 
years. 

No ellipsoidal variations have been observed in RS 
Oph, which is generally assumed to indicate that it 
does not fill its Roche lobe. If this is the case, then 
the accretion disk must b e fed by some mechanism 
other than RLO. Work by Ppdsiadlowski & Mo hamedl 
(2007) explores one such alternative. Those authors pre- 
sented a scenario dubbed "wind Roche lobe overflow", 
where the stellar wind of the giant fills the Roche lobe, 
even though its photosphere (determined from mass- 
luminosity-radius relationships) does not. The wind is 
gravitationally focused towards the LI point, where it is 
then captured by the potential well of the white dwarf. 
At this point, the material can form an accretion disk, as 
in RLO systems. This mode of mass transfer can be up 
to 100 times more efficient than Bondi-Hoyle accretion. 
If such mass transfer is occurring in RS Oph, it can feed 
a disk, and the white dwarf, at a sufficiently high rate to 
power the observed outburst frequency without the red 
giant filling its Roche lobe. 

In this paper, we present new deep, quiescent-state 
X-ray observations of RS Oph, carried out with Chan- 
dra and XMM-Newton 537 and 744 days after the 2006 
outburst, respectively. These observations shed light on 
the origins of the X-ray emission in the system, and re- 
veal new details about the accretion of material onto the 
white dwarf component. In Section [2] we describe the 
observations and the data reduction techniques used to 
obtain lightcurves and spectra for analysis. We present 
our timing analysis in Section [3] and our spectral analy- 
sis and best fitting models in Section [4] We discuss the 
derived quiescence properties of RS Oph in the context 
of earlier X-ray observations of the system in Section [5j 
In Section [6j we discuss our findings and their implica- 
tions for understanding the physics of accretion in the 
RS Oph system. Finally, a summary of our work is given 
in Section 

2. OBSERVATIONS AND DATA REDUCTION 

The first of our observations of RS Oph was carried 
out on 2007 August 4 (537 days after the 2006 outburst 
began) with the ACIS-S camera onboard the Chandra 
X-ray observatory, for a total exposure time of 91 ks 
(see Table 1). This exposure was triggered as part of 
a target of opportunity (ToO) program to image jet ac- 
tivity in symbiotic stars. Faint extended str ucture is in 
fact p resent in the data, and is presented in iLuna et al.l 
(2009). The exposure was performed in FAINT mode, 
with the source positioned on the ACIS-S3 chip (one of 
two back-illuminated CCDs) which has the highest sen- 
sitivity to soft photons. The data as downloaded from 
the archive had been reduced with the most recent cal- 
ibration files available, so no further reprocessing of the 
data was necessary. We extracted timing and spectral 
products using the Chandra Interactive Analysis of Ob- 
servations (CIAO) software, version 4.1.1. 
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Table 1 

New observations of RS Oph 



Date Observatory Instrument Net t exp (s) Days from outburst a Flux (erg s 1 cm 2 ) b 

2007- 08-04 06:11:26 Chandra ACIS-S 90993 537 5.2 x 10" 13 

2008- 02-26 16:57:59 XMM-Newton EPIC-pn 20780 744 2.1 X 10~ 13 

EPIC-MOS1 28168 

EPIC-MOS2 28461 

a We set t = 2,453,779.33 after Narumi et al. (2006). 

b Observed flux, not corrected for absorption. Fluxes are in the energy range 0.3—10.0 keV 



Source lightcurves and spectra were extracted from a 
6" circle centered on the source using the CIAO task 
dmextract. Background data were extracted from an 
annulus also centered on the source, with inner and outer 
radii of 10" and 20" respectively. There was a total of 
6468 good events in the source region, with 264 back- 
ground counts, or 4% of the observed counts. Spec- 
tral response files were generated using the CIAO tasks 
mkacisrmf (response matrix) and mkarf (effective area). 
In order to facilitate the use of \ 2 statistics, we grouped 
the spectra to give a minimum of thirty counts per bin. 
Using PIMMS, we verified that our instrumental setup 
resulted in a pile up rate of less than 2%. 

We observed RS Oph for a second time on 2008 Febru- 
ary 26 (day 744 after outburst), this time with the 
XMM-Newton observatory. All instruments onboard 
were turned on during the observation. The European 
Photon Imaging Cameras (EPIC) MOS and pn instru- 
ments were all operated in full window mode, which ex- 
poses the entire field of view of the camera. We reduced 
these data using the Science Analysis Software (SAS) v. 
11.0. First, we reran the pipeline routines in SAS in or- 
der to apply the most up to date calibration to the data. 
The observation, which had a total duration of 35 ks, 
was affected by two periods of high background - one at 
the start of the observation, and a shorter interval near 
the end. We screened these periods of high background 
from the event file by defining a new good time interval 
(GTI) file, which excludes time periods where the EPIC- 
pn (EPIC-MOS) 10-12 keV count rate is >0.4 (0.35) cts 
s _1 . This reduced the total exposure time for each instru- 
ment (see Table 1), and also resulted in a low number of 
counts above 2 keV. In order to improve the statistics for 
spectral fitting at high energies, we filtered the good time 
interval data by grade, retaining both single and double 
events i.e. PATTERN < =4, in the EPIC-pn datfl For 
the MOS data, we retained events with grades <= 12. 
We masked bad pixels and CCD edge defects using the 
most conservative screening criterion, FLAG==0, as rec- 
ommended for spectral analysis. 

Source lightcurves and spectra were extracted from a 
30" circular region for all three EPIC instruments. Back- 
ground events in the two MOS datasets were extracted 
from an annular region of inner radius 35" and outer ra- 
dius 100" . For the pn data, the background events were 
extracted from two rectangular regions at the same chip 

7 In general, it is recommended that only events with grade 
(i.e. singles only) are retained for spectral analysis of EPIC-pn 
data. Using both singles and doubles can in principle increase 
the pile up rate, and reduce the energy resolution of the resulting 
spectrum. However RS Oph does not have a high enough count 
rate for this effect to be a concern. 



y coordinate, chosen to maximize the background area 
while extracting events with similar distance from the 
read out register, as recommended in the XMM-Newton 
EPIC calibration documentation. A total of 2047 good 
counts were extracted in the source region in the EPIC- 
pn data, with 74 of these (4%) due to the background. 
The number of counts extracted from the MOS data 
was much lower - 630 (636) from MOS1 (MOS2), with 
119 (109) background counts. Spectral response files for 
all the EPIC spectra were created using the SAS tasks 
rmf gen and arf gen. Again, we binned the spectra with 
a minimum of thirty counts per bin. With PIMMS, we 
verified that the pile up rate in this case was less than 
1%. 

The two reflection grating spectrograph (RGS) instru- 
ments were operated in standard spectroscopy mode, but 
no usable data was obtained due to the low count rate 
and we do not discuss them further. The optical moni- 
tor (OM), operating in imaging mode, obtained images 
of the source in the V, U and UVM2 filters, with most 
exposures lasting 800s (the default for the standard mo- 
saic pattern). The complete dataset is comprised of 5 
images in each of the V and U filters, and 11 in the 
UVM2 filter. Note that the OM can also be operated 
in fast mode, in which the detector is continuously read 
out every 0.5s. Use of this mode would be extremely use- 
ful in any future observations of RS Oph, since it pro- 
vides much higher time resolution than standard imaging 
mode. We analyzed the pipeline reduced images, which 
include photometry of the brightest sources in each im- 
age. The pipeline carries out the following photometric 
procedure: source counts are extracted from a 16.5 arc- 
sec radius circle, and background counts from a 18.3-20.7 
arcsec annulus. These are converted to count rates by di- 
viding by the exposure time. Next, these count rates are 
corrected for bad pixels and photon coincidence losses, 
which are well understood for the OM. Finally, back- 
ground counts are subtracted from the source annulus, 
based on the corrected count rates. Magnitudes are then 
calculated using this corrected count rate and the zero- 
point defined for each filter in the XMM documentation^ 

3. TIMING ANALYSIS 

The Chandra lightcurves are presented in Figure 1, 
binned into 4000 s segments. The most striking feature 
of these lightcurves is the rise in count rate during the 
first half of the observation, which increases from 0.05 
to 0.08 c s _1 in the 0.3-10 keV energy range. We ex- 
tracted X-ray lightcurves from each observation across 
the entire energy response of the instrument, and in three 
smaller energy bands: 0.3-1, 1-2 and 2-10 keV. Looking 

8 http:/ /xmm.esac.esa.int/sas/current/howtousesas.shtml 
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Figure 1. ACIS-S lightcurves in various energy ranges, binned in 
4000s intervals. The vertical dashed line indicates the time used 
to split the observation into two spectra (see text). 
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Figure 2. XMM-Newton EPIC-pn lightcurve, binned in 1000s 
intervals. 



at the lightcurve behavior in the smaller energy bands, 
it is clear that the count rate increase occurs primarily 
at intermediate energies, between 1 and 2 keV (see Fig. 
1, green data points). The XMM-Newton lightcurve is 
shown in Figure 2. There are relatively more counts in 
the soft band compared to the medium and hard bands 
than in the Chandra observation, although some of this 
effect can be explained by the different energy response 
of the EPIC-pn and ACIS-S instruments. However, that 
the 0.3-10 keV count rates are similar shows that the 
system must be intrinsically fainter during the XMM- 
Newton observation, since XMM-Newton has a much 
larger collecting area than Chandra. We do not observe 
any significant count rate changes over the course of the 
XMM-Newton observation, although we note that expo- 
sure time is shorter than the timescales over which such 
changes occurred in the Chandra observation. 

The OM lightcurves in the U, V and UVM2 filters 
are shown in Figure 3, with no correction for interstellar 
reddening. The U — V color is con s istent with the val- 
ues reported bv iHenden fc Munaril (pOOl . The UVM2 
data show low- level variability from exposure to expo- 
sure. Following ISokoloski et al.l (|200lD . we compared the 
standard deviation of the pipeline measured magnitudes 
with that expected if the observed changes in magnitude 
were due only to the errors (Poisson plus systematic) in 
the individual measurements. This is the same as test- 
ing the hypothesis that the observed lightcurve is consis- 
tent with a source of constant flux using the % 2 statis- 
tic. The expected standard deviation, s exp , is found from 
the average of the observed errors on each measurement. 
These are calculated by the XMM SAS pipeline, based 
on the Poisson errors of the observed source and back- 
ground counts. Taking the average value is appropriate 
here since the exposure time of each observation was the 
same. 

We find that the observed standard deviation of the 
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Figure 3. Optical Monitor V, U and UVM2 lightcurves. Magni- 
tudes are in the standard Vega system. Note that the uncertainties 
on each point are smaller than the size of the symbols. 



UVM2 brightness (0.077 mags) is ~2 times larger than 
that expected from the measured errors alone (0.037 
mags), a firm detection of variability. This could be 
due to accretion induced flickering, which tends to be 
enhanced in bluer filters where the disk light dominates 
over the emission from the cooler companion. If we work 
under the assumption that the observed fluctuations are 
due to flickering, and hence accretion, we can use the 
amplitude of the flickering to place constraints on the ac- 
cretion ra te two years a fter the outburst. We follow the 
method of lBruchl(p92T ). who assumed that the flickering 
source was 100% modulated and that all other sources 
of light were constant. If the amplitude of the flicker- 
ing is Am and the mean observed magnitude over the 
course of the observations is rh, then the magnitude of 
the constant source, m c is given by 



0.5Ar 



(2) 



Assuming that the faintest data point in the lightcurve is 
due only to the flux of the constant source, F c , and that 
the brightest data point is due to the constant source plus 
the flickering source (F c + Ff), then it is straightforward 
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Table 2 

Accretion rate lower limits from UV flickering 

1000 



E(B-V) A 23 io m F f (erg s" 1 cm" 2 ) M (M 0> r/y yr" 1 ) 

0.4 3.51 11.05 1.5 x 10 -11 2.3 x lO" 10 

0.5 4.38 10.18 3.4 x 10" 11 5.3 x 10" 10 

0.6 5.26 9.30 7.4 x 10" 11 1.1 x 10" 9 

0.7 6.14 8.43 1.7 x 10~ 10 2.6 x 10~ 9 
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Figure 4. Chandra ACIS-S spectrum with model CA-va. The 
unbinned data have been smoothed for the purposes of plotting to 
give a minimum signal to noise ratio of 3 in each bin. Dotted lines 
indicate the individual components of the model, and the ratio of 
the model to the data is shown in the lower panel. 



prised of a number of spectral lines. The most obvious 
of these can be identified as the unresolved He-like line 
triplets of Mg XI, Si XIII and S XV at -1.3, 1.8 and 2.4 
keV, respectively. We also observe an emission feature 
at —6.5 keV, which is likely the unresolved blend of some 
combination of the neutral iron fluorescence line at 6.4 
keV, the He- like Fe triplet at 6.7 keV, and the Lyman a 
line of H-like Fe at 6.9 keV. 

Finally, there is a large, sharp decrease in continuum 
flux at energies higher than 7 keV. This absorption fea- 
ture is due to the 7.1 keV K edge of iron, and is an 
unmistakable signature of strong absorption of the hard 
X-ray continuum. We can estimate the column density 
of Fe by comparing the count rate on either side of the 
edge (equivalent to flux since the instrumental effective 
area is approximately constant over this energy range) 
and using the following relationship: 

N(Fe) = -ln{F out /F in )/a FeK (4) 

where F out is the flux after the edge, Fi n the flux be- 
fore the edge (excluding the emission line feature), and 
tJFeK the cross sec t ion o f the Fe K edge. Following 
ILeahv fe Creightonl (|1993lh we adopt a apeK of 3.6 x 
10~ 20 cm 2 . In the Chandra spectrum, the flux drops 
by a factor of —5 at the Fe K edge, giving N(Fe) = 
—4 x 10 19 cm" 2 . This can be converted to a hydro- 
gen column, N(H), by multiplying by the ratio of the 
numbers of H and Fe at oms in a solar abu n dance gas. 
Using the abundances of iGrevesse fe Sauvall (|1998l ). the 
observed Fe edge implies a column density N(H) — 1.4 x 
10 24 cm~ 2 . High column density absorption has been 
observed in other symb iotic stars, including RT Cru 
dLuna fe Sokoloskil 12007ft and T CrB, another recurrent 
nova (|Luna et al.ll2008f ). Since such a high column den- 
sity absorber removes nearly all flux below — 2 keV (in- 
deed, the spectra of both RT Cru and T CrB have almost 
no flux at low energies), we should expect any model of 
the emission in RS Oph to be comprised of multiple com- 
ponents subject to different absorption in order to simul- 
taneously account for the large soft X-ray flux observed 
below 1 keV, and the deep Fe K edge. 



a The value quoted is the lower limit to the accretion rate, since no 
bolometric correction is made to the derived accretion luminosity. 
The error on the flux measurement (and therefore the accretion 
rate) due to uncertainties in the extinction law towards RS Oph 
is ~60%. The mass accretion rate values assume a distance of 1.6 
kpc. 



to show that 
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We obtain both Am (0.23 mags) and m (14.56 mags) 
from the lightcurve. The resulting ratio of Ff /F c is 0.24. 

In order to estimate the actual flux of the flickering 
source, the observed magnitude must first be corrected 
for the effects of interstellar extinction. A number of 
different E(B-V) estimates are given in the literature, 
ranging from 0.4 to 0.73. This spread in the estimate of 
E(B-V) is the dominant source of uncertainty in calcu- 
lating the luminosity of the flickering source, since this 
corresponds to a difference in the magnitude correction 
of 2.63, or a change in the measured flux of a factor of 
— 10. Deciding which value to use depends somewhat on 
where the UV light is believed to originate. Some intrin- 
sic absorption of the accretion disk is likely, given that 
the white dwarf is embedded within the wind of the red 
giant. Therefore, we calculated the flickering source flux 
for a range of reddening values. First, we obtained the 
extinction coeffici ent in the UVM2 filte r, A2310, using the 
reddening law of lCardelli et all (|1989f ). for E(B-V) val- 
ues of 0.4, 0.5, 0.6 and 0.7. These values are presented 
in Table 2, along with mo, the reddening corrected mean 
magnitude for each value of E(B-V). While the spread in 
E(B-V) is the largest source of uncertainty, an important 
source of error for a given E(B-V) value is the potential 
deviation from the assumed Car delli et al. extinc tion 
law. From the data presented in iFitzpatr ick (199^), we 
estimate that this effect contributes a 0.5 mag error into 
the extinction, and thus into the corrected mean mag- 
nitude. This corresponds to an uncertainty in the de- 
rived flux of —68%. Next, we converted the ifiQ values 
into count rates, and then to fluxes, using the conversion 
factors outlined in the Optical Monitor calibration doc- 
umentation, and the effective width of the UVM2 filter, 
400 A. Finally, we obtained the flux of the flickering com- 
ponent by multiplying by the ratio found using Equation 
3. These fluxes are presented in Table 2, and are in the 
range 1.5 x 10~ n to 1.7 x 10 -10 erg s _1 cm -2 , depend- 
ing on the assumed value of E(B-V). 

4. SPECTRAL ANALYSIS 
4.1. Chandra 

The high quality of the Chandra data, presented in 
Figure 4, reveals a number of notable spectral features 
that can guide us in developing a model of the emission. 
We observe soft emission, peaking around 1 keV, com- 
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We develope d a model fo r the Chandra data using 
Xspec vl2.0.5 JA rnaudlfl996l) and the binned spectrum 
in the energy range 0.3-10.0 keV. All quoted uncertain- 
ties are the 90% confidence intervals obtained using the 
error command in Xspec. If during the execution of 
this command a better fit was obtained (based on the % 2 
statistic), the error analysis was restarted using the new 
best fit parameters. This process was repeated, iterating 
on the parameter values, until a stable model was found 
and the error command completed without finding a new 
fit. Throughout this section, we quote distance depen- 
dent quantities in units of D = 1.6 kpc (jHicllming et al.l 
1986), which was originally derived from HI 21 cm radio 
data. This distance is the one most commonly used in 
the literature, and is in agreement with the value 1.4+2 \ 
kpc derived by IBarrv et al.l (|2008j ) from a survey of all 
available distance estimate methods. We will discuss the 
implications of using larger distances in Section 6.3. 

We used t he ap ec models in Xspec, developed by 
iSmith et al.l ([200 Is), which model both continuum and 
line emission from thermal plasmas. We found no accept- 
able fits for solar abundance plasma models, even with 
up to three temperature components subject to different 
absorbing columns. This is not surprising, since there is 
clear evidence of non-solar abundances in the RS Oph 
system from both outburst and quiescence observations. 
Based on analysis of high resolution grati ng spectra of 
RS Op h obtained during the 2006 outburst . INelson et al.l 
( 2008) and INess et al.l (|2009[ ) both concluded that the 
X-ray emitting materi al in RS Oph i s enha nced in ni- 
trogen. Furthermore, iPavlenko et al.l (|2008l ) concluded 
from near-infrared spectroscopy of CO and CN molecu- 
lar features that the wind from the red giant is depleted 
in carbon and enhanced in nitrogen, possibly due to ei- 
ther CNO processing in the hydrogen-burning shell of the 
red giant or contamination from earlier nova outbursts. 

Therefore, in order to explore the effects of non-solar 
abundances on our model fits, we make use of the values 
determined by Ness et al. (2009) from the June 2006 
Chandra LETG observation of RS Oph. This observa- 
tion occurred after the end of the supersoft source stage, 
when the X-rays likely originated in the shocked ejecta, 
and therefore are the most appropriate to compare with 
the quiescence data. This set of abundances, taken from 
Table 7 in Ness et al., is presented in the first column of 
Table 3. The values presen t ed ar e relative to the solar 
values of lGrevesse fc Sauvall ((19980 , and we assume solar 
values for elements not listed in Table 3. It is important 
to note that the Ness et al. abundances were determined 
using an emission measure distribution method that as- 
sumed a solar abundance of oxygen, which may or may 
not be appropriate for emission from the nova ejecta. We 
implemented these abundances in Xspec using the vapec 
model, which allows the abundance of each element to be 
varied individually. 

We find no acceptable fit (x 2 jv — 2.9 for v = 101) for 
a two plasma model, where each component is subject 
to differing absorption. Adding a third absorbed plasma 
component improves the fit, particularly at the soft end 
of the spectrum, with x 2 l v — 2-1 for v = 98, although it 
is still not statistically acceptable. The soft flux arises in 
two components with kT = 0.2 and 0.6 keV, subject to 
absorbing columns of 4.7 and 1.0 x 10 21 cm -2 , respec- 



Table 3 

Elemental abundances used in model components 
Value relative to solar a 



Element 


N09 


vapec 


O 


1.00 b 


0.7±0.1 


N 


7.54^ 2 


7.54 


Ne 


1.36 ± 0.16 


1.36 


Mg 


1.91 ± 0.70 


1.5 ± 0.3 


Si 


1.62 ± 1.00 


1.62 


Fe 


0.43 ± 0.10 


0.6 ± 0.1 



vmcf low 



average 



7+0.2 
-0.1 
7.54 
1.36 
1 6+ - 4 
1.62 



0.7Z 



0.6 1 



7+0.2 
-0.1 

7.54 

1.36 

1 55+ ' 4 

1.62 
q +0.2 
-0.1 



0.7Z 



0.6 1 



a Solar abundances are taken from Grevesse and Sauval (1998). 
b These abundances were derived using the Emission Measure Dis- 
tribution (EMD) method described in Ness at al. (2009), which 
assumes a solar abundance of oxygen. 



tively. These values of N(H) are about a factor of ^2 
higher and lower, respectively, than the i nterstellar value 
of 2.4 ± 0.6 x 10 21 cm -2 reported by [Hiellming et al.l 
(1986) from HI 21cm radio observations. The flux above 
2 keV originates in a third component, whose tempera- 
ture in this model is 69 keV. However, the fit to the data 
(particularly above 5 keV) is extremely poor, and the 
model has a much steeper slope at these energies than 
is observed in the data. Also, the column density of the 
absorber attenuating this component is extremely small 
(3.7 x 10 19 cm -2 ), which is completely at odds with the 
observed Fe K edge at 7.1 keV. The low column den- 
sity of this absorber is required by the model in order 
to reproduce the continuum slope between 2 and 5 keV 
- a complete absorber with a column density of ~10 24 
cm -2 removes all 0.3-2 keV, and results in a rising con- 
tinuum energies greater than 2 keV, which is clearly not 
observed in the data. This discrepancy, considered in 
tandem with best fit N(H) values which differ signifi- 
cantly from the interstellar value towards RS Oph, made 
us search for alternative models which can account for 
the spectral properties in a more consistent way. 

A possible answer to the question of the physical ori- 
gin of the X-ray emission of RS Oph may lie in X-ray 
studies of other symbiotic stars. The fits to the X-ray 
spectra of RT Cru and T CrB both require the pres- 
ence of an absorber that only partially covers the X-ray 
emitting source to explain their 2-10 keV emission. Al- 
though difficult to disentangle in RS Oph (where there 
is clearly an additional, less absorbed soft component) 
we also tried to fit the continuum emission with a model 
comprised of a completely absorbed soft component, and 
a partially covered, absorbed hard component. This has 
the effect of allowing soft flux from the hard component 
to contribute to the spectrum, even though the column 
in front of the hard X-ray emitting material is very high. 
We used the pcf abs model in Xspec for this case, which 
is characterized by an absorbing column N(H) pc , and a 
covering fraction (C. F. < 1) that quantifies the degree 
to which the absorber covers the source. 

The fit statistic for this model (which we implemented 
as wabs*vapec+pcf abs*vapec in Xspec) is much better 
than the model with a completely covered hard compo- 
nent, with x 2 l v — 1-5 for v = 119. The temperature 
of the soft component is 0.35^q;o4 keV subject to an ab- 
sorbing column N(H) tot = 2.4t°; 2 x 10 21 cm -2 , which 
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is consistent with the ISM contribution. The hard com- 
ponent has a temperature of 3.1^Q'g keV. The absorber 
has a column density N(H) pc = 1.3 ± 0.2 X 10 24 cm~ 2 , 
and covers 98 ± 1 % of the source. There is a narrow, 3a 
residual feature at approximately 6.4 keV, which likely 
indicates the presence of the Fe K alpha fluorescence line. 
This line is produced via fluorescence after ionization of a 
K shell electron. Given the high optical depth of the Fe K 
edge, the K a line should be bright. We therefore added 
a gaussian to the model at 6.4 keV, and kept the width 
consistent with instrumental broadening only. The best 
fit energy of the line is 6.4 ± 0.04 keV, and the equivalent 
width is 320±2?o 

eV. All other quantities related to the 
harder plasma and absorber remain the same within the 
uncertainties as the model with no gaussian. 

Examining the residuals of the model fit to the data, 
we see that the model predicts more flux than is ob- 
served at energies below 0.5 keV. All of this excess 
flux arises in the partially covered component of the 
model. However, even this component should be sub- 
ject to interstellar absorption along the line of sight. 
Since the absorption affecting the soft component is 
consistent with being due only to the ISM, we al- 
tered the model slightly so that all components are af- 
fected by a complete absorber, representing the ISM, 
and the hard component further attenuated by the par- 
tial covering absorber. In Xspec, this was implemented 
as wabs (vapec+pcf abs (vapec+gauss) ) . Although this 
modification makes physical sense, it does not formally 
improve the fit, giving a x 2 l v — 1-5 for v = 100. In fact, 
the model now underpredicts the flux at energies less 
than 0.6 keV. As a final refinement, we added an addi- 
tional plasma absorbed only by the total absorber, imple- 
mented as wabs (vapec+vapec+pcf abs (vapec+gauss) ) 
in Xspec. This results in an improved reduced \ 2 value 
of 1.3 for v = 98. The first two plasma components have 
temperatures of 0.14tg;°2 and 0.35±°;q 2 keV, absorbed 

by a column with N(H) tot = 2.lt°A x 1q21 cm" 2 , which 
is still consistent with previous estimates of the ISM con- 
tribution towards RS Oph. The third plasma component 
has kT = 2.5 to\l keV, and the partial covering absorber 

has N(H) pc = l.l±o.2 x 1024 crrr2 ' with a 98 ± 1 % 
covering fraction. 

We note that we have binned the data in order to utilize 
the x 2 statistic in developing our model. Such binning 
unavoidably leads to a loss of spectral information. In or- 
der to check if our parameter values were in anyway influ- 
enced by our choice of binning, we refit o ur model us ing 
the unbinned data and the Cash statistic (Cash 1979) as 
implemented in XSpec. We find no difference in the best 
fit model parameters using the different statistic, and the 
range of uncertainties is also similar, showing that our 
model fitting is not influenced in any strong way by the 
choice of binning. However, since the fit using the Cash 
statistic makes use of all the available spectral informa- 
tion, we show the parameters and uncertainties derived 
using that statistic in Table 4, as model CA. We plot the 
data with this model in Figure 4, and have smoothed the 
data for plotting purposes to give a minimum signal to 
noise ratio in each bin of 3. 

4.1.1. A cooling flow origin for the hard component 



The X-ray spectra of many accreting white dwarfs, 
including the hard X-ray detected symbiot ic stars, can 
be w ell described by a cooling flow model (Muk aiet al.l 
2003). In this model, which was originally developed for 
the gas in galaxy clusters, hot material in the post shock 
region of the accretion flow cools isobarically an d settles 
onto the white dwarf surface (sec Fabian (1994) for a re- 
view of the cooling flow model). The model assumes no 
additional sources of heating or cooling once the accreted 
material is in the boundary layer. While this is generally 
the case for non-magnetic white dwarfs (for magnetic sys- 
tems, cyclotron cooling can be important), it may not in 
fact be true for RS Oph, as we will discuss in Section 6. 
The model is particularly useful for studying accretion in 
CVs, since it provides an estimate of the accretion rate 
onto the white dwarf directly from model fitting. To 
determine whether the RS Oph data can be described 
by this model, we replaced the partially absorbed vapec 
component in model A with vmcf low, the variable abun- 
dance implementation of the cooling flow model in Xspec. 
We used the Ness et al. (2009) abundances, and allowed 
the maximum and minimum temperatures of the cooling 
flow to vary freely. 

The resulting model has an equally good fit to the data 
as the vapec model (x 2 '/v = 1.29, v — 95). We checked 
the fit result using the unbinned data and the Cash statis- 
tic, and again found no differences in best fit parame- 
ter values within the uncertainties. We present the best 
fit parameters derived using the Cash statistic for this 
model (model CB) in Table 4. The parameters for the 
soft component (i.e the two vapec components) agree 
with those of model CA within the uncertainties. The 
parameters of the total and partial covering absorbers 
are also the same within the uncertainties. We find min- 
imum and maximum temperatures for the cooling flow 
of <0.97 keV and 5. 91^2 keV, respectively. The upper 
limit to the minimum temperature indicates that cooling 
within the flow is efficient. The maximum temperature 
in the cooling flow is surprisingly low, and we will discuss 
this result further in Section 6.2. Finally, the normaliza- 
tion gives an estimate of the cooling flow accretion rate, 
in this case 2.4± 2 '2 x 10~ 9 (D/1.6 kpc) 2 M yr~\ where 
D is the distance to the source in kpc. 

4.1.2. Effects of different abundance? 

We explored whether allowing the abundances of our 
two best fit models to vary would improve the fit. This 
was motivated primarily by the fact that the abundance 
estimates of Ness et al. (2009) have an associated un- 
certainty. Determining abundances from CCD resolution 
spectrum is difficult, since few spectral features are avail- 
able for use in model fitting. In this case, we allowed the 
abundances of Ne, Mg, Si and Fe to vary freely, since 
these spectral features are clearly visible in the data. We 
also explored the effect of changing the O abundance, 
since a clear residual at 0.6 keV (due to O VIII Ly a) 
is observed in the fit using the Ness abundances (which 
were determined assuming solar oxygen, see Section 4.1). 
For each of these elements, we identified the model com- 
ponent contributing the majority of the flux, and allowed 
the elemental abundance parameter in that component 
to vary freely. The associated abundance parameters in 
the other model components were tied to this value. By 
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Table 4 

Model fits to the Chandra spectra 







Complete 


observation 




First 


half 


Second half 


Parameter 


CA 


CB 


CA-va 


CB-va 


CIA 


C1B 


C2A 


C2B 



N(H) totai (10 21 cm- 2 ) 
kTi (keV) 
normi (10 — 4 ) 
kT 2 (keV) 
norrri2 (lO -4 ) 
N(H) pc a (10 22 cm- 2 ) 
C. F. b 
kT 3 (keV) 
norm3 (10 -3 ) 
kT mln (keV) 
kT maI (keV) 

M/D'i 6kpc (10-9 M yr-l) 

E gaua3 (keV) 

norm 9 „„ ss (lO -6 ) 



1.9 



+0.4 
0.4 



n 1 ^+ - 06 

u - io -0.03 

o.4±g:| 
0-3718:°; 

2-61S1 
117±|S 
0.98+ 0.01 

n c+0.5 
z -°-0.3 

7l| 



1 

0.16 



+0.3 
-0.4 
+0.08 

0.04 

o.3l°;| 

o-37i8:8? 

n q+0.5 

120l 4 ° 
0.98±0.01 



<1.8 

5 4+ 2 ' 4 
J - -l.i 

2.6+ 2 - 7 



2.1+ 0.5 
0.25 4 
2.^ 



h0.04 
-0.02 
1-1.2 
-0.7 

1.3+ 05 



-o.t 

8128 

0.98± 0.01 



h0.6 
-0.5 
h0.02 
-0.03 
hi. 3 
-0.7 
n co+0.09 

u-oj-o.os 

1 1+0.5 



2.1! 
0.25^ 
2.1* 



106! 



108! 



a li 

0.97+0.02 



2 6+°' 7 
z -°-0.4 

c+4 

°-2 



6 3O+0-03 
"• oy -0.04 
1+9 
-5 



6.39 
11 



1.9 
+0.03 
0.04 
11 



6.39 
9.3 



0.03 
-0.04 
+6.9 
4.3 



<1.2 
5.9l?j 

i-6i : 9 9 

6.39+0.03 

^+8.8 
-3.8 



9.4^ 



0.26! 



2 

-+0.02 
J -0.05 

2-2l?:? 

0.43+0.13 



0.9 



+ 1.8 



-0.5 
4+30 
-19 

0.98+ 01 



84! 



2.4 



0.02 
+0.8 
-0.5 
4+5 



^+0.03 
J -0.07 
1+3.3 
-1.1 
,+0.13 
-0.25 
1+1.6 
-0.7 

0.9718:8^ 



0.26! 
2.2^ 

0.43j 
0.9l 
86"* 



+0.02 
0.03 



2 

0.25 

1.9+0.3 
0.60+0.04 

1 5+°' 4 

11311 
0.97+ ,' 01 



0.25+0.03 
1.9+0.3 
60+ 04 



1.5 



+0.3 



2.6 
5 



6.39 

1Z ' Z -5.0 



1.6 
3.7 
3.6 
6 

11 



1.7 
1.5 
+ 14.4 

1.7 
+46.4 
1.9 

6.39 
n 

5 



O.03 
+0.7 
0.5 
+6 
3 



-0.5 
d+60 
3 -40 

07+0.02 

U - M '_0.04 



118! 



4.7 



6.39 

+9.2 

4.7 



<3.7 

4 2+ 4 ' 8 
^• z -1.7 
+ 110.9 

3.2 

6.39 

+8.0 
3.8 



3.8 



4.9 1 



Fc 
O 

Mg 

C-statistic (f) 



0.43 
1.0 
1.91 



0.43 
1.0 
1.91 



0.6+ 0.1 
0.7+0.1 
1.5+ 0.3 



0.6 
0.7 
1.6 



1.27 (96) 
773 (652) 



1.29 (95) 
775 (651) 



1.03 (93) 
742 (649) 



o.i 

+0.2 
0.1 
+0.4 
-0.3 



0.6 
0.7 
1.55 



0.6 
0.7 
1.55 



0.6 
0.7 
1.55 



0.6 
0.7 
1.55 



1.04 (92) 
743 (648) 



0.90 (53) 0.93 (52) 0.84 (75) 1.04 (92) 
608 (654) 610 (653) 619 (654) 621 (653) 



Note. — "A" models are of the form wabs (vapec+vapec+pcf abs (vapec+gauss) ) , i.e. have a single temperature plasma as the hard 
component. "B" models are of the form wabs (vapec+vapec+pcf abs (vmcf low+gauss) ) i.e. have a cooling flow as the hard component. 
Models CA and CB assume the abundances of Ness et al. 2009, as outlined in Section 4. In models CA-va and CB-va, the abundances of 
Fe , O and Mg were allowed to vary freely (see Section 4.1.2). The separate fits to the first (CI models) and second half (C2 models) of the 
observation use the average abundance values found in models CA-va and CB-va. All quoted uncertainties arc 90% confidence intervals 
obtained using the error command in Xspec. Parameters with no uncertainties were not allowed to vary, as discussed in the body of the 
text. 

a The column density of the partial covering absorber. 
a The covering fraction of the partial covering absorber. 



doing this, all spectral components varied in tandem, but 
the model fitting was guided by the component in which 
the spectral feature of interest was strongest. 

We initially varied each elemental abundance individ- 
ually, starting from either model CA or CB, and used 
the F-test to determine whether changing the abundance 
parameter resulted in a statistically significant improve- 
ment in the reduced \ 2 value. We find that allowing Fc, 
O or Mg to vary freely results in an improved fit that is 
significant at a level greater than 99.9%. The change in 
the Fe abundance is the most significant in terms of the 
F-test result. Freeing the Ne abundance improves the fit 
marginally, but the F-test null hypothesis can only be 
rejected at the 75% level, so the improvement is prob- 
ably not significant. Varying Si results in no change to 
the fit statistic. We find similar improvements in the fit 
to both the single temperature and cooling flow models 
developed above. 

Next, we refit the models allowing the O, Mg and Fe 
abundances to all vary freely. The resulting models have 
smaller \ 2 l v values and residuals when compared to the 
Ness et al. (2009) abundances. Once again, fitting the 
unbinned data using the Cash statistic did not change the 
best fit parameters within the uncertainties determined 
using the \ 2 statistic. We present the resulting abun- 
dances obtained for each model (and their average val- 
ues) in Table 3, and the model parameters (model CA-va 
with vapec hard component; model CB-va with vmcf low 
hard component, where "-va" indicates that abundances 
were varied) obtained using the C-statistic method, in 



Table 4. The best fit values agree within uncertainties 
between the two models. We actually determine a sub- 
solar abundance of O, in contrast with the assumption 
of solar O in the DEM method of Ness et al.. However, 
we note that our best fit value (0.7l^i) is m good agree- 
ment with the O abundance determined by Ness et al. 
from fitting apec models to the Day 26 grating data, a 
method that estimates absolute abundance corrections 
by comparing line strengths with the continuum. We 
find a higher Fc and lower Mg abundance than Ness et 
al. (0.618 ? an d l-55l8'3i respectively), although within 
the uncertainties they are in agreement. The best fit 
temperatures and normalizations of the two soft plasmas 
are different than in models CA and CB (see Table 4), 
but the hard component parameters remain in agreement 
within the uncertainties. 

4.1.3. Count rate and spectral changes during the Chandra 

observation 

In order to explore the physical origin of the count rate 
variations observed in the Chandra lightcurve, we ex- 
tracted two separate spectra from the observation — one 
from the first half, where the count rate was increasing, 
and one from the second half where it was approximately 
constant. The time dividing the observation is marked 
in Fig 1 by the dashed vertical line. We used the best fit 
models found for the total spectrum as the starting point 
for fitting the split spectra. Since looking at smaller time 
periods reduces the signal to noise of the spectrum, we 
have fixed the values of a number of parameters in or- 
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Figure 5. Comparison of the spectra obtained from the first half 
(black data points) and second half (red data points) of the Chandra 
observation. We also show the associated models CIA and C2A. 
The data have been smoothed as in Figure 4. The ratio of the 
model to the data for each spectrum is shown in the lower panel. 

der to reduce the number of degrees of freedom. The 
N(H) of the total absorber (N(H) tot ) is consistent with 
the ISM contribution in all models, and therefore we set 
N(TL) total = 2 x 10 21 cm~ 2 . We set the Fe, O and Mg 
abundances to the average values from models C and D 
i.e. 0.6, 0.7 and 1.55 times the solar values. Finally, we 
set E gauss to 6.39 keV. Reducing the degrees of freedom 
in this reasonable way makes a stable fit more likely, and 
allows us to explore which parameters are responsible for 
the change in count rate. The resulting best fit parame- 
ters and associated uncertainties are listed in the second 
half of Table 4. The spectra and models CIA and C2A 
are presented in Figure 5, with the first half of the ob- 
servation in black, and the second in red. Overall, the 
two spectra are quite similar — within the uncertainties in 
each bin, they are the same below 0.8 keV. They are also 
similar above 2 keV, although the Fe Ka feature at 6.4 
keV is brighter in the first half than in the second. The 
main differences between the two spectra occur between 
0.8 and 2.0 keV, as the lightcurves indicated. 

However, we cannot definitively identify the origin of 
this change from our best fit models — within the for- 
mal statistical errors there is no difference between the 
two halves of the spectrum. Comparing the two halves 
using the x 2_s tatistic, the parameters of the soft plas- 
mas agree within uncertainties, and all changes appear 
to be in the hard component (both in the column den- 
sity of the absorber, and in the normalization of the hard 
plasma). However, checking the fits using the unbinned 
data and the C-statistic, we find instead that the change 
in 1-2 keV count rate is due to a change of temperature 
and normalization in the second soft plasma component, 
with no detectable difference in the hard plasma. 

4.2. XMM-Newton 

The high background level during parts of the XMM- 
Newton observation resulted in a net exposure time of 
only ~21 ks for the pn camera, and ~28 ks for the MOS 
cameras. Thus there are only a small number of counts 
above 2 keV, which increases the challenge of determining 
the best fit model to the data. We fit our model to the 
spectra obtained with all three cameras simultaneously, 



and we used the models developed for the Chandra data 
in Table 4 as the starting parameter values. As we did for 
the split spectra discussed in Section 4.1.3, we fixed the 
N(H) tot (the column density of the total absorber) to 2 x 
10 21 cm -2 , and the energy of the gaussian component to 
6.39 keV to decrease the number of degrees of freedom. 
Again, since the total absorber is consistent with the ISM 
contribution, which should not evolve with time, we are 
justified in constraining this parameter. We used two 
sets of abundances - those of Ness et al. (2009), and 
the updated and averaged Fe, O and Mg abundances we 
found from the Chandra observation, to see if there were 
differences in the fit parameters, as we found with the 
Chandra data. 

The resulting best fits to the XMM-Newton data are 
presented in Table 6 (models XA, XB, XA-n and XB-n), 
and the data and model XA-va are presented in Figure 
6. Here, the "n" models use the average new abundances 
found from the fits to the Chandra spectrum in this pa- 
per, like the fits to the partial spectra in Table 4. The 
emission below 2 keV is still well described by a two 
temperature plasma model with kT values between 0.14 
and 0.62 keV. In contrast with the Chandra analysis, 
there is no clear difference in the soft plasma best fit pa- 
rameters between the two sets of abundances. The tem- 
perature values agree with those found in the Chandra 
observation within the uncertainties. However, the nor- 
malization of these components (and thus the emission 
measure of the plasmas) has decreased. While a hard 
component is still clearly seen in the data, the parame- 
ters are poorly constrained due to the small number of 
counts. The best fit values for the normalization (single 
temperature) or accretion rate (cooling flow) of the hard 
component are much lower than for the Chandra obser- 
vation, although the formal confidence intervals are large 
and do not completely rule out higher values. This large 
uncertainty is driven by degeneracies in the N(H) pc , cov- 
ering fraction and plasma normalization space. With so 
few hard counts, there are few features in the data that 
can eliminate these degeneracies, as we were able to do in 
the Chandra data. However, the lower observed flux and 
smaller hard component normalizations suggest that the 
accretion rate is lower than in the Chandra observation. 

4.3. Constraints on an optically thick boundary layer 

Several authors have shown that at high accretion 
rates, the boundary layers in CVs should become op- 
tically thick, thermalize, and emit like a blackbody with 
effective temperature of a few tens of eV (depending on 

M); see e.g. iPatterson fe Ravmondl ()1985[ ). This black- 
body component would be subject to the same absorp- 
tion as the hard component. Therefore, in order to place 
constraints on the luminosity and temperature of an op- 
tically thick boundary layer, we added a blackbody com- 
ponent to our best fit Chandra and XMM-Newton mod- 
els. The luminosity of this optically thick portion of the 
boundary layer is given by 

Lbl = 47r R 2 WD o- T^ff Ibl (5) 

where Jbl is the fraction of the white dwarf sur- 
face covered by the optically thick boundary layer, ex- 
pected to be <10% ( see e.g. iPopham fe Naravanl 119951 : 
iPiro fe Bildstenll2004ft . In our model we assume Rw d = 
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Table 5 

Model fits to the XMM-Newton spectrum 



Parameter 



Model XA 



Model XB 



Model XA-n 



Model XB-n 



N(H) totai (10 
kTi (keV) 
normi (10 — 4 ) 
kT 2 (keV) 
norm2 (10 — 4 ) 
N(H) pc a (10 22 
C. F. b 
kT 3 (keV) 
norm3 (10~ 3 ) 
kT mi „ (keV) 
kT max (keV) 

E 9 ouss (keV) 
norm flatlss 



cm 2 ) 



2 

,+0.03 
-0.02 
,+0.5 
-0.4 

0.31 ± 0.01 

7+0-1 
-0.2 
,+36 
-58 
,+0.03 
-0.04 
=+9 
3 -3 
,+ 1.0 
-0.4 



0.12: 

.31 A 
1.7j 

nr 

0.96j 
6j 
0.6"* 



0.12" 1 
l.F 1 



2 

,+0.02 
-0.03 
,+1.0 
-0.5 

0.31 ± 0.02 

7+0-1 
-0.2 

,+54 
-44 



1-7Z 
109^ 



0.95 



+0.03 
-0.03 



0.13 
0.9 

0.30 
1.9 
124 

0.95 
10 
0.6 



M yr" 1 ) 



6.39 



2.3 



+7.2 



<11 
14+47 

1+ 21 
Uil -0.07 

6.39 

,+6.9 
-1.7 



2 

+0.04 
0.03 
+0.5 
0.3 
+0.02 
0.01 
+0.2 
0.3 
-73 
-45 
+0.03 
0.06 
+34 
4 

+ 1.0 
0.3 



0.13 
0.8 

0.29 
1.9 
119 

0.95 



2 

+0.04 
0.03 
+0.8 
0.3 
+0.02 
0.01 
+0.2 
-0.3 
+ 72 
43 
+0.03 
0.05 



2.2"* 



2.1 



6.39 

+6.7 
2.0 



< 4.4 
>12 

+0.22 
0.07 

6.39 

n q+5.8 
z - z -1.9 



0.09 



Fe 


0.43 


0.43 


0.6 


0.6 


O 


1.0 


1.0 


0.7 


0.7 


Mg 


1.91 


1.91 


1.55 


1.55 


X z /v O) 


1.61 (85) 


1.64 (84) 


1.42 (85) 


1.44 (84) 


C-statistic (f) 


2732 (5811) 


2732 (5810) 


2734 (5811) 


2732 (5810) 



Note. — The form of these models is the same as that outlined in Table 4. Again, parameters with no uncertainties were fixed as 
outlined in the text. The "-n" models denote that the abundances determined in this paper were used. 
a The column density of the partial covering absorber. 



The covering fraction of the partial covering absorber. 




" 1 , or an accretion rate of 5.7 x 10 
-l 



Figure 6. XMM-Newton EPIC-pn spectrum with model XA-n. 
The pn data are in black, and the MOS1 and MOS2 data in red 
and green, respectively. We show the ratio of the model to the data 
in the lower panel. The data have been smoothed as in Figure 4. 



3 x 10 8 cm, /bl = 0.1, leave T e // as the free parame- 
ter, and tie the normalization constant of the blackbody 
component (the luminosity) to T e f f via the above equa- 
tion, assuming a distance of 1.6 kpc. The blackbody 
component is absorbed by both the interstellar medium 
and the partial covering absorber in our model. We find 
that the maximum YT e ff allowed in our model is 0.034 
keV (~395,000 K) for the Chandra model, and 0.031 
keV (-360,000 K) for XMM-Newton. For temperatures 
higher than this, a significant excess begins to emerge 
over the data at energies <0.5 keV. This result depends 
somewhat on the partial covering fraction of the complex 
absorber — the higher the covering fraction, the higher al- 
lowable temperature. 
The maximum allowable blackbody temperature value 



of 0.034 keV for the Chandra observation corresponds 
to an intrinsic luminosity of 1.6 x 10 35 (D/1.6 kpc) 2 
erg s _1 . If this was the only source of emission from 
the boundary layer, and assuming Mwd — 1-3 M Q , we 
would obtain a value for M of 8.5 x 10~ 9 (-D/1.6 kpc) 2 
M Q yr -1 . For the XMM-Newton data, the maximum 
temperature implies an accretion luminosity of 1.1 x 10 35 
(-D/1.6 kpc) 2 erg s" 
(D/1.6 kpc) 2 M yr~ 

The derived blackbody temperature limits indicate 
that most of the flux from any optically thick boundary 
layer would emerge in the UV. Therefore, as a consis- 
tency check, we predicted the flux in the OM UVM2 fil- 
ter using one of our best fit EPIC-pn model (model XA) 
plus the maximum allowable blackbody (31 eV, 51.1 x 
10 35 erg s _1 ). We did this by creating a one channel 
spectrum file using the average observed UVM2 count 
rate, and included this data point in our XSpec model by 
using response files downloaded from the XMM-Newton 
calibration database. We find that the UV flux predicted 
by the model is a factor of a few hundred lower than is 
actually observed. Thus, adding the blackbody remains 
consistent with the OM data. 

5. COMPARISON WITH PREVIOUS X-RAY 
OBSERVATIONS 

We compared the best fit models and fluxes from the 
current data to those obtained by ROSAT and ASCA 
before the 2006 outburst. The observed count rates of 
the old observations are low, and therefore the statistics 
are poor. However, by comparing our best fit models 
and observed count rates at different energy bands, we 
can investigate the long term behavior of the different 
emission components observed in the system. In order 
to estimate equivalent ROSAT and ASCA count rates, 
we folded our best fit Chandra and XMM-Newton mod- 
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els through the appropriate instrumental responses in 
XSpec. These count rates were then compared to the 
previously published data. 

The ROSAT observations are useful for looking at the 
soft component of the emission. In our best fit models, 
this arises in warm plasma with temperatures ranging 
from 0.15 to 0.4 keV. The first ROSAT observation, car- 
ried out on 1991 March 2, had an observed count rate 
in the energy range 0.2-2.4 keV of 3.5 x 10 -3 c/s. This 
had increased to 1.2 ± 0.2 x 10~ 2 c/s one year later, on 
1992 March 4. If the source emission was described by 
our best fit Chandra models, all of which result in simi- 
lar fluxes, then the resulting ROSAT 0.15-4.0 keV count 
rate would have been 2.1 x 10~ 2 c/s, a factor of two 
higher than the brightest ROSAT observation. Our best 
fit XMM-Newton model is in better agreement with the 
data, differing by only 25% from the observed count rate 
in the 1993 observations. However, in 1991 the system 
was much fainter than during our XMM-Newton obser- 
vation. 

Our best fit Chandra and XMM-Newton models are 
both consistent with the 1997 ASCA/GIS data in the 2- 
10 keV energy range, primarily due to the poor counting 
statistics. The observed GIS count rate in the 2.0-10.0 
keV range is 1.2 ± 0.8 x 10~ 3 c/s, compared with 1.4 
x 10~ 3 c/s (model CA, Chandra) and 3.2 x 10~ 4 c/s 
(model XA, XMM-Newton). The Chandra models over- 
estimate the flux in the 0.7-2.0 keV range, as observed 
with the ASCA/STS instrument. However, the XMM- 
Newton models are in much better agreement with the 
data. 

6. INTERPRETATION OF THE DATA AND MODEL FITS 
6.1. Evidence of high column intrinsic absorption 

All of our best fit models are comprised of a soft and 
hard component. The hard component is equally well 
described by either a single temperature plasma, or a 
cooling flow, plus a 6.4 keV neutral Fe emission line ob- 
scured by a high column density intrinsic absorber which 
covers most, but not all, of the hard X-ray emitting re- 
gion. Extremely high values for the column density of 
this complex absorber are implied by the presence of a 
deep Fe K edge at 7.1 keV, and the values of this param- 
eter in our best fit models are between 0.9 and 1.2 x 10 24 
cm~ 2 . In all models, the covering fraction is >95%. Sim- 
ilar high column intrinsic absorbers have been required to 
fit the X-ray s pectra of other symb iotic stars; T CrB, RT 
Cru, CH Cyg (|Kennea et al.ll2009h . and perhaps most in- 
terestingly the jet prod ucing system MWC 5 60, which is 
viewed almost pole-on (Stutc & Sahai 2009). Since only 
the hard component is affected by this absorption, the 
a bsorbing mate r ial m ust be close to the white dwarf. 

Waldcr et al.1 (|2008D carried out hydrodynamical mod- 
eling of an accretion-outburst cycle in RS Oph, assum- 
ing a mass loss rate from the red giant of 10 -7 M Q yr _1 . 
They find that a thick disk forms with a strong spiral dis- 
tortion, and that the accretion rate onto the white dwarf 
is 10% of the mass loss rate from the giant. The densities 
in the vicinity of the white dwarf are as high as 10 -12 ' 5 
g cm -3 , which depending on the line of sight through 
the material can result in column densities as high as a 
few 10 23 cm -2 , consistent with the values implied by our 
observations. Therefore, it is possible that the accretion 



disk plays some role in absorbing the X-rays from the 
boundary layer. 

The existence of a high column density absorber in the 
system offers a way to reconcile a high accretion rate with 
a low observed X-ray luminosity. Correcting only for 
interstellar absorption, the luminosity of the hard com- 
ponent in the 0.3-10 keV range is 1.1 x 10 32 (D/1.6 
kpc) 2 erg s _1 . However, when we also correct for the 
presence of the complex absorber, the luminosity in the 
same range becomes 1.8 x 10 33 (D/1.6 kpc) 2 erg s _1 — 
a factor of 10 increase. Our limits on an optically thick 
component of the boundary layer show that the unab- 
sorbed luminosity of the system could be as high as 1.6 
x 10 35 (D/1.6 kpc) 2 erg s _1 . Therefore, the high degree 
of intrinsic absorption in the system can easily reduce the 
emergent X-ray flux by a factor of 1000. This is one fac- 
tor that explains the so-called "missing boundary layer" 
problem in RS Oph. 

6.2. A low temperature accretion component 

Given the high mass of the white dwarf in RS Oph, 
the low maximum temperature of the cooling flow sug- 
gests that some other cooling mechanism besides thermal 
bremsstrahlung cooling is at work. Models of optically 
thin, radial accretion onto white dwarfs predict that the 
maximum temperature of the post shock gas is set by 
the freefall velocity, and therefore the mass, of the white 
dwarf as follows: 

3^rn H GM WD . . 

J- shock = 57-5 (0) 

where /i is the mean molecular weight (assumed here to 
be 0.63), ran the mass of one hydrogen atom, and Mwd 
and Rwd the mass and radius of the white dwarf re- 
spectively. For accretion through a disk, the velocity of 
the pre-shock gas is given by the innermost Keplerian 
orbit, and is a factor of V2 less than in the free fall case. 
Therefore the resulting temperature, which varies as v 2 , 
is a factor of two lower. Assuming AfwD = 1.3 Mq and 
i?WD = 3 x 10 8 cm for RS Oph, the expected T s h oc k 
is ~63 keV. Cooling flows with high maximum tempera- 
tures are seen in several other symbiotic stars believed to 
harbor white dwarfs as massive a s the one in RS Oph (e.g. 
RT C ru (jLuna fe Sokoloskil2007f) and T CrB (jLuna et all 
2008)). However, the maximum temperature of the flow 
in RS Oph derived from our model fits (< 6 keV where 
the statistics are good) is much lower than in those sys- 
tems. We note that the mass accretion rates implied by 
the cooling flow normalization in our models are of the 
same order as those in T CrB and RT Cru, i.e a few 10 -9 
M Q yr" 1 . Why is the maximum temperature so much 
lower? 

One possible explanation is that cooling by Compton 
scattering is important in the post shock accretion flow 
in RS Oph. In this scenario, there is an intense pho- 
ton field that Compton scatters off, and cools, the hot 
post-shock electrons. The overall plasma temperature 
will drop if the energy losses due to inverse Compton 
scattering d o minat e over those due to free-free emission. 
iFrank et all (|2002j ) (hereafter FKR02) give the following 
expression for the ratio of the Compton cooling time to 
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the free-free cooling time (their Equation 6.56): 
tcompt _ 7.5 x lCT 5 iV e 



1/2 



(7) 



where N e and T e are the electron density and temper- 
ature of the emitting plasma and U ra d is the radiation 
energy density of the scattering light source. Assuming 
that the seed photons arise from some source on the sur- 
face of the white dwarf, U ra d can be approximated as 
follows. 

Urad — ~. ; (8) 

^Rw D fradC 

Here, L ra d is the luminosity of the radiation source pro- 
viding the scattering photons, f ra d is the fraction of the 
white dwarf surface over which this light is emitted, and 
c is the speed of light. The post shock density is deter- 
mined by the mass continuity equation and the strong 
shock jump conditions so that 



M 



C 



■VK 



(9) 



where m e is the electron mass, f acc is the fraction of the 
white dwarf surface over which mass is accreted (note 
that this is not necessarily the same as fsti the fraction 
of the white dwarf surface taken up by an optically thick 
boundary layer, in Equation 7), and vk is the innermost 
Keplerian velocity. Inserting equations 9 and 10 into Eq. 
8, expressing T e using Equation 7 in the Keplerian case, 
and noting that vk — {GMwd/ Rwd)®' 5 , we find that 
the ratio of the cooling times can be expressed as follows: 



tCompt 



- 2 - 3 L rad,33 M -9 R 9 °' 5 M 1.3 fradf, 



-1.5 . 



r 1 . 

' acc 



(10) 



Here, we have expressed the luminosity of the radiation 
source, L ra( j,33, in units of 10 33 erg s~ x . M_g is the mass 
accretion rate onto the white dwarf in units of 10 -9 M Q 
yr _1 , i?g is the radius of the white dwarf in units of 10 9 
cm, and Mi. 3 is the white dwarf mass in units of 1.3 Mq. 

In order for Compton cooling to be important in RS 
Oph, we require that tcompt/tff < E If we assume that 
the post shock flow is optically thin, and that the low 
temperature is due to Compton cooling, then the cool- 
ing flow model will reflect the correct accretion rate, but 
have a maximum temperature representative of the re- 
gion in the flow where the density increases to the point 
where free-free emission once again dominates the cool- 
ing. Therefore, we can take M from the cooling flow 
models developed to describe the data. The other pa- 
rameters, f ra d and /acc, are unknown. / acc depends on 
the structure of the inner disk and boundary layer, which 
is model dependent. An estimate can be made by calcu- 
lating the ratio of the scale heigh t of the boundary layer 
to the radius of the white dwarf. iPatterson fc Raymond! 
(1985) give an expression for the boundary layer scale 
height, which we have expressed in terms of our white 
dwarf parameters and accretion rates as follows: 



H = 5.1 X 10 7 M 1 °3 1 Mrg' 5 cm, 



(11) 



where Mi. 3 and M_g have the same meaning as in Equa- 
tion 11. We find that on day 537, the average M obtained 



from the cooling flow component in our best fit models 
is ~2 x 10 -9 M Q yr^ 1 , resulting in a corresponding / acc 
of 0.1 (we concentrate on this epoch since the day 744 
XMM-Newton observation does not allow us to place con- 
straints on the accretion rate) . The fraction of the white 
dwarf surface which emits the scattering photon field, 
frad, depends on the physical origin of the radiation. 

One source could be residual photospheric emission 
from the white dwarf, in which case f ra d — 1- During 
a nova outburst, the outer layers of the white dwarf are 
heated by the nuclear burning taking place on the sur- 
face. Once burning ceases, the white dwarf must radiate 
away this excess heat. In work exploring the evolution of 
a clas sical nova throu gh a complete accretion-outburst 
cycle, iPrialnikl (|1986l ) found that shortly after nuclear 
burning ends, the white dwarf experiences a dramatic 
drop in luminosity by a factor of 10 as the nuclear ashes 
cool, and then a longer phase of luminosity decline as the 
white dwarf loses the excess heat in its outer layers (see 
their Figure 8). During that phase, which lasts about 
10 years, L oc t~ 1A4 , where t is the time since the on- 
set of the outburst. For this radiation source to cool the 
post shock plasma in RS Oph, we find that L ra d > 8.4 
x 10 34 erg s~ 4 . This implies an effective temperature of 
the white dwarf photosphere of ^190, 000 K. Taking the 
post-burning luminosity of the white dwarf given in Pri- 
alnik (1986) of 10 36 erg s -1 , then at 1.5 years after the 
outburst we expect the residual luminosity to be of order 
10 35 erg s _1 . This would provide the required radiation 
field to cool the post shock accretion flow. We note that 
such a source would be undetectable in our data, even if 
it is only absorbed by the interstellar medium. 

Alternatively, if the boundary layer is partially opti- 
cally thick, then this additional component could pro- 
vide the necessary radiation field to cool the optically 
thin portion of the post shock plasma. If this is the case, 
then frad would be determined by the structure of the 
optically thick part of the boundary layer. As we dis- 
cussed in Section 4.3, an optically thick boundary layer 
is expected to cover ~10% of the white dwarf surface. 
Again taking an average M in the optically thin part of 
the boundary layer of 2 x 10~ 9 M Q yr _1 , then an / acc 
of 0.1 requires L ra d > 8.4 x 10 33 erg s^ 1 for cooling 
due to Compton scattering to be dominant. The limits 
we placed on the existence of an optically thick bound- 
ary layer show that such a source of radiation could be 
present without being detected directly in the observed 
spectrum of the system. 

We acknowledge that invoking Compton scattering as 
the explanation for the low cooling flow maximum tem- 
perature requires a source of radiation which we have 
not directly observed in the data. Therefore we also con- 
sider a scenario that does not require the presence of 
such a radiation field. It may be that RS Oph falls into 
a range of accretion rates where the transition from an 
opticall y thin to an optically thick boundary layer takes 
place. IPopham fc Naravanl (|1995f ) explore the proper- 
ties of optically thick boundary layers for white dwarfs 
with masses in the range 0.6-1.0 M Q , particularly their 
dependence on accretion rate, white dwarf mass and ro- 
tation rate. They find that the transition from optically 
thin to thick occurs at higher accretion rates for more 
massive white dwarfs (see their Figure 8, upper panel). 
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They note that their studies indicate a transition region 
between the two states, where the boundary layer is op- 
tically thin, but has not yet become hot. For a one solar 
mass white dwarf, this would occur for M between 3 x 
10 -10 and 10 -9 M yr _1 . However, the transition should 
occur at higher M for a more massive white dwarf. Un- 
fortunately, the authors do not elaborate on this regime 
further - for now we just note that this could offer an 
explanation for the lack of a very high temperature com- 
ponent in the RS Oph spectrum. 

Any of these scenarios can potentially explain the dif- 
ference between RS Oph and T CrB, another recurrent 
nova and symbiotic star. The longer interval between 
nova outbursts in T CrB than in RS Oph implies a lower 
accretion rate in T CrB. If this is the case, then these 
tw o systems may strad d le the transition zone discussed 
in iPopham fc Narayanl (|1995f ). thus explaining the dif- 
ferent cooling flow maximum temperatures. A slightly 
lower accretion rate may also mean that T CrB has ei- 
ther no optically thick boundary layer component, or if 
one is present, it may not be luminous enough to result 
in tcompt/tf f < 1. If the white dwarf photosphere is the 
main source of Compton scattering photons in RS Oph, 
then the longer time since the last outburst in T CrB (60 
years vs. 2 years) and the lower expected photospheric 
luminosity may result in less efficient Compton cooling 
of the post shock flow, and the higher observed maxi- 
mum temperature. A potential test of this hypothesis 
is to look for a low cooling flow maximum temperature 
within the first few years of the next outburst of T CrB. 

6.3. The mass accretion rate onto the white dwarf 

At each epoch, the spectral models we have developed 
allow us to explore the range of possible accretion rates 
onto the white dwarf. As we discussed in Section 1 , exist- 
ing estimates of the accretion rate in RS Oph using X-ray 
observations are too low to account for the observed out - 
burst frequency. Using the results of lYaron et al.l (|2005h . 
the triggering mass required to power an outburst every 
~20 years is 3.5 x 10" 6 M Q for a 1.2 M Q WD, and 3.47 
x 10~ 7 M Q for a 1.4 M WD. The associated accretion 
rates are therefore 1.8 x 1CT 7 and 1.7 x 1CT 8 M© yr -1 , 
respectively. For each dataset, we have found two equally 
well fitting models. In the first, the hard component is 
a single temperature plasma, and in the second it is a 
cooling flow. 

We can estimate the accretion rate from the single tem- 
perature model by obtaining the unabsorbed X-ray flux 
of the hard component, converting this to a luminosity, 
and inserting into Equation 1 (again assuming a 1.3 M 
white dwarf with radius 3 x 10 8 cm). Note that the lu- 
minosity depends on the distance to RS Oph, and there 
is a degree of uncertainty in this quantity. Throughout 
this paper, we have assumed a distance of 1.6 kpc (see 
Section 4.1). However, larger dist ance estimates exist in 
the literature. iRupen et all (|2008l ) proposed D = 2.45 ± 
0.4 kpc, based on the size of the radio outflow spatially 
resolved with the Very L ong Basel i ne Arr ay (VLBA) dur- 
ing the 2006 outburst. Schacfei] (|2010D favors an even 
larger distance, ~4 kpc, based on an average of a num- 
ber of different methods. However, at this distance the 
nebular expansion observed in both the radio and the 
optical with the Hubble Space telescope (jRibeiro et al.l 



[200^1 imply expansion velocities of >8000 kms 1 , which 
appears to be too large for a nova. For this reason, we 
take the larger value of Rupen et al. as the upper limit 
to the distance. 

We obtain the unabsorbed flux of the single temper- 
ature hard plasma in XSpec by setting the normaliza- 
tion of the softer plasma components, and all absorption 
components, to zero. An associated uncertainty is de- 
rived by looking at how the flux changes within the 90% 
confidence interval range for both the plasma tempera- 
ture and normalization. For the Chandra observation, 
we find a range of unabsorbed 0.3-10 keV fluxes for this 
component of 0.6-2.0 x 10 -11 erg s _1 cm" 2 for both sets 
of abundances (i.e. models CA and CA-va in Table 4). 
Assuming a distance of 1.6 kpc, these flux measurements 
imply accretion rates of 1.1-3.1 x 10~ 10 (D/1.6 kpc) 2 
Mq yr -1 . In the XMM-Newton observation the flux val- 
ues are a factor of ~10 lower (0.3-3.2 x 10~ 12 erg s _1 
cm -2 ), and hence so is the implied accretion rate (0.6- 
4.5 x 10" 11 (D/1.6 kpc) 2 M yr" 1 ). Placing RS Oph 
at 2.45 kpc (the Rupen et al. distance) increases the X- 
ray luminosity, and hence accretion rate, by a factor of 
2.3. Thus, the maximum accretion rate is 7.1 x 10~ 10 
M Q yr -1 1.5 years after outburst, and 10 -10 M yr -1 2 
years after outburst. 

The cooling flow model provides a direct estimate of 
the mass accretion rate onto the white dwarf, if we as- 
sume that the boundary layer is completely optically thin 
(if an optically thick region exists the accretion rate can 
be higher, as we'll discuss shortly). The values found 
for the Chandra observation are in the range 0.7-5.3 x 
10~ 9 Mq yr -1 — factors of 3-50 times larger than those 
inferred from the single temperature plasma model. Ac- 
counting for the uncertainty in the distance, the accre- 
tion rate implied by the cooling flow fits can be as large 
as 1.2 x 10~ 8 using the radio distance of Rupen et al. 
This value is within a factor of 2 of the rate required to 
power an outburst in the system every ~20 years, if the 
white dwarf is close to 1.4 M Q . The estimate obtained for 
the XMM-Newton observation is an order of magnitude 
lower, although with much larger uncertainties given the 
poorer statistics of that data. In any case, it seems clear 
from the variations in X-ray brightness between the two 
epochs that the accretion rate is variable in quiescence 
by up to a factor of 10. 

If the boundary layer is partially optically thick, which 
may be the case in RS Oph as we discussed above, then 
the data are consistent with an even higher mass accre- 
tion rate. As we showed in Section 4.3, a boundary layer 
with luminosity up to 1.6 x 10 35 (15/1.6 kpc) 2 erg s _1 is 
consistent with the Chandra data, assuming our best fit 
model is the correct one. This implies an accretion rate 
of 8.5 x 10" 9 (D/1.6 kpc) 2 M Q yr" 1 through an optically 
thick component of the boundary layer. Adding this to 
the M estimate from the cooling flow, and placing the 
system at 2.45 kpc results in an accretion rate can power 
an outburst every 20 years. 

We can use the UV flickering to place a further con- 
straint on the accretion rate 2 years after the outburst. 
Using the values of the flux of the flickering source (Ff) 

presented in Table 2, a lower limit to M (since we have 
not performed any bolometric correction and therefore 
have not accounted for all the flux) can be obtained us- 
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ing the following equation 



M > ^ D * RwD Ff 



(12) 



GMwd 

Again assuming a distance to RS Oph of 1.6 kpc, a white 
dwarf mass M\yd of 1.3 M©, and a white dwarf radius 
Rwd of 3 x 10 8 cm, we can express M in a more conve- 
nient form 

M > 2.3 x lO" 10 ( — — — ] I — ] 

M\YD 



1.3 M 



1.6 kpc J V3xl0 8 cm 



1.5 x 10 -11 erg s 1 cm 2 



) Moyr- 1 . 



(13) 



We present the values of M at this epoch for our range 
of assumed reddening in Table 2. We favor an E(B-V) 
for the UV light of 0.7 , since this value was derived 
from earlier UV observations obtained with IUE during 
both outburst and quiescence ()Sniiderslll987|) . This E(B- 
V) value corresponds to an N(H) = 3.7 x 10 21 cm~ 2 , 
assuming a conversion factor N(H ) = 5.3 x 10 21 E(B- 
V) cm -2 (|Predehl fe Schmittlll995f) . This value is higher 
than the ISM absorption, most likely due to intrinsic 
attenuation of the UV flux by the red giant wind. For 
this value of the reddening, the corresponding lower limit 
to the mass accretion rate is 2.6 x 1CP 9 M yr _1 . This 
is much higher than the accretion rate estimates given by 
Orio (1993) and Mukai (2008), and is much closer the the 
expected accretion rate given the frequency of observed 
outburst, as discussed in lYaron et al.l (|2005l ). 

Ultimately, we have shown that if the hard X-ray emis- 
sion originates in a cooling flow, then the implied accre- 
tion rate is very close to the theoretical rate given by nova 
theory if the mass of the white dwarf is within ^10% of 
Chandrasekhar limit. If the boundary layer is partially 
optically thick, then the data are consistent with the high 
accretion rate if RS Oph is at the distance of Rupen et 
al. (2.45 kpc). Since the model fits suggest that the ac- 
cretion rate is variable in quiescence, then it is possible 
that at other times the accretion rate is higher. Thus, 
when averaged over the inter-outburst time period, the 
required M is reached. If the cooling flow model is incor- 
rect, then we must seek some other mechanism to supply 
the white dwarf with mass, since the rates inferred from 
the single temperature plasma models are much too low 
(even at the maximum distance) to result in the observed 
outburst frequency. This would most likely be an opti- 
cally thick boundary layer, which is hidden from our view 
by the ISM and the intrinsic absorber in the system. 

6.4. Emission from the nova shell 

The soft component is generally described by a two 
temperature plasma with non-solar abundances. The 
temperatures and relative fluxes vary slightly between 
our models, but the presence of a less absorbed, low tem- 
perature plasma is evident from a direct examination of 
the data (since such soft flux is at odds with the large Fe 
K edge observed at 7.1 keV). This component is observed 
to fade between the Chandra and XMM-Newton observa- 
tions, but the temperatures agree within the uncertain- 
ties. Therefore, it is likely that this component arises 



in the shocked ejecta and red giant wind created during 
the early stages of the outburst. The derived emission 
measures (obtained from the normalization constants of 
the thermal plasma components, assuming a distance of 
1.6 kpc) are in the range 10 54 -10 55 cm -3 — a factor of 
1000-100 00 lower than observe d immediately after the 
outburst (Sokolosk i et al.ll2006l ). The cooling time for a 
shocked optically thin plasma can be calculated by di- 
viding the total thermal energy by the energy lost per 
unit volume 



^cool — 



2 ^e^max 

^A(T) 



(14) 



where n is the electron density, T max the maximum post 
shock temperature and A(T) the radiative cooling func- 
tion at a given T. Over the range of temperatures where 
the plasma emits in the X-rays, the cooling function is ap- 
proximately co nstant and has a value of ~10~ 23 erg s _1 
cm 3 (see e.g. iRavmond fe Smith! 119771 ). Taking T max 
= 10 s K and an estimate of n e in the envelope of the 
red gi ant obtained 2-3 yrs after t he 1985 outburst of 10 7 
cm~ 3 (lAnupama fc P rabhu 1989), we find a cooling time 
of the order of 10 yrs (this n e estimate is of the same order 
of magnit ude as many o t her sy mbiotic stars as demon- 
strated in lLuna fc Costal (|2005l i). Thus, it seems reason- 
able that the shocked material would still be emitting 
at the times of the new observations. The two temper- 
ature fit is most likely an approximation to the actual 
temperature distribution of the plasma, which should be 
continuous. Our best fit XMM-Newton model can ac- 
count for the ROSAT emission observed 7 years after 
outburst, which suggests that the shell emits for a large 
fraction of the inter-outburst time. The soft count rate 
variation observed between the first and second ROSAT 
observations is larger than that seen over the course of 
the Chandra observation (factor of 3 for ROSAT versus 
2 for Chandra), but they may both be due to the same 
effect i.e. some combination of changing absorption of 
the hard component. 

The new observations also allow us to examine the long 
term evolution of the system post outburst. The X-ray 
flux has continued to decrease since the end of the su- 
persoft phase, and is a factor of ^2 lower in the energy 
range 0.3-10 keV in the day 537 observation than in the 
XMM-Newton observation obtained 239 days after the 
outburst . By the time of the XMM-Newton pointing on 
day 744, the flux is a factor of ^6 lower than on day 239. 
This change of flux occ urs primarily in the soft compo- 
nent. INess et al.l (|2009D compared the 1.1-1.8 keV flux 
in all available Chandra, XMM-Newton and Swift ob- 
servations obtained after the 2006 outburst (this small 
energy band was chosen to avoid the supersoft emission 
which appeared several weeks into the outburst), They 
found that after June 2006, the flux decayed as t~ 8 / 3 , 
where t is the time since outburst. To see if the new 
data are consistent with this rate of decline, we com- 
pared the observed count rates in the same energy range 
to the rates predicted by the < -8 / 3 decay, using the 1.1- 
1.8 keV EPIC-pn count rates on days 205 and 239 of the 
outburst as our initial values. The observed flux is a fac- 
tor of 5 and 2 higher at 90% confidence level than the 
values predicted for days 537 and 744, respectively. This 
is likely due to the reestablishment of accretion — in our 
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models, the harder plasma component is the dominant 
source of flux above 1.4 keV. 

6.5. Comparison of UV results with optical studies 

Recently, iZamanov et al.l (|2010f ) published a study of 
optical flickering in RS Oph, using optical photometry 
obtained in July 2008 and July 2009. Based on the de- 
gree of flickering in each of their program filters (U, B, 
V or I), they calculate the flux of the flickering compo- 
nent at multiple wavelengths. They then deredden these 
values assuming E(B-V) = 0.73. Finally, they determine 
the best blackbody fit to the derived flickering spectrum 
at the two observation epochs. They find that the source 
of the optical flickering has a temperature of around 9500 
K in both observations, but that the flickering source is 
significantly brighter in 2009 than in 2008, with a bolo- 
metric luminosity of ~150 L© (vs. 50 L© in 2008). These 
values correspond to accretion rates (assuming the same 
white dwarf parameters as we did in Section 4.3) of 3.2 
and 1.5 x 10 -8 M© yr _1 assuming a distance of 1.6 kpc. 

Although our UV observations were carried out before 
both of these observing epochs, it is useful to compare 
our measured flux with their best fit models of the optical 
flickering source. To facilitate this comparison, we note 
that Zamanov et al. define the flickering flux, Ff as F - 
F rnin , whereas we define F c = F max - F min , resulting in a 
difference in flux values of a factor of ~2. Assuming the 
same degree of reddening as Zamanov et al. (i.e. E(B-V) 
= 0.73), we obtain a flickering source flux density from 
our OM UVM2 (2310 A) photometry of 5.3 (±3.2) x 
10~ 13 erg s^ 1 cm~ 2 A -1 . This is higher than the value 
predicted by their best fit to the optical data in July 2008, 
and consistent with the July 2009 fit (see their Figure 5). 
The temperature found for the optical flickering source is 
much too low to be associated with the boundary layer, 
and is more likely to be a signature of the accretion disk. 

7. CONCLUSIONS 

We have analyzed X-ray and UV data from two obser- 
vations of RS Oph carried out 537 and 744 days after the 
2006 outburst. The Chandra observation reveals that the 
X-ray emission is variable on short timescales, with the 
changes most apparent at energies between 1 and 2 keV. 
A comparison of the two datasets also reveals long term 
variability, with the 0.3-10 keV flux reduced by a factor 
of 3 in the XMM-Newton observation. When compared 
with the long term brightness evolution observed with 
Swift, we find that RS Oph is brighter than the trend 
observed in the Swift data would predict, indicating the 
existence of a new source of flux in the system, most likely 
renewed accretion. The X-ray spectra are well described 
by two components - a soft optically thin plasma compo- 
nent with temperatures 0.15-0.65 keV, absorbed only by 
the interstellar medium, and a harder, highly absorbed 
(N(H) ~ 10 24 cm -2 , 97-99% covering fraction) compo- 
nent, which is equally well described by a single temper- 
ature plasma or a cooling flow. The soft component is 
observed to decrease in brightness between the two ob- 
servations, although the temperature remains constant 
within our measured uncertainties. We propose that this 
emission arises in the shocked nova ejecta, which is still 
cooling from its peak post-outburst temperature, and can 
remain an X-ray emitter for up to a decade after the out- 



burst. The hard component arises in the accretion disk 
boundary layer. 

Our best fit models can also be used to describe the 
emission observed with ROSAT 6 and 7 years after the 
1985 outburst, and with ASCA 12 years after the out- 
burst. The current data have brighter soft emission than 
was observed with either satellite in the quiescence pe- 
riod before the 2006 outburst. This is most likely because 
in those observations the nova shell emission, which dom- 
inates below 1.4 keV, was much fainter due to the larger 
amount of time passed since outburst. The models pre- 
sented here for both Chandra and XMM-Newton can be 
consistent with the hard emission observed with ASCA, 
primarily due to the poor statistics of that dataset. 
ROSAT was not particularly sensitive to the accretion 
component, which is most obvious above 2 keV. How- 
ever, the differences in count rate observed with ROSAT 
could also be explained by a change in the hard compo- 
nent, as we saw during the Chandra observation on day 
537. 

The temperature associated with the hard component 
is only 2.61q4 keV (for the single temperature plasma) or 
5.9l 3 ;g keV (for the cooling flow model). This is surpris- 
ingly low for accretion onto such a massive white dwarf, 
and strikingly different from the maximum temperatures 
found in similar systems like T CrB. We can conceive 
of two explanations for this low temperature value. The 
first is that cooling via Compton scattering dominates 
the accretion flow, in the boundary layer, lowering the 
maximum temperature of the emission. The scattering 
radiation field could originate in either the white dwarf 
photosphere, which was recently heated by the nuclear 
burning associated with the nova outburst, or from some 
optically thick region in the boundary layer. The ab- 
sence of such a source in T CrB (either due to a lower 
accretion rate or the longer time since the last nova out- 
burst) would explain the difference in observed maximum 
temperature. The second possibility is that M is such 
that the system lies in the transition region between a 
completely optically thick and completely optically thin 
boundary layer. In this case T CrB and RS Oph straddle 
the transition accretion rate, and therefore have different 
spectral characteristics. 

Our model fits to the X-ray spectra obtained on days 
537 and 744, and the presence of UV flickering during the 
day 744 observation with XMM-Newton, can be used to 
estimate a range of reasonable accretion rates onto the 
white dwarf. In the X-rays, the high quality Chandra 
data provide the most robust estimates of M. Uncer- 
tainties in the accretion rate come from the fact that 
two physical models describe the data equally well, and 
also from the range of distances proposed for the system. 
Spectral changes over the course of the Chandra obser- 
vation on day 537, and the lower observed 2-10 keV flux 
on day 744, strongly suggest that the accretion rate onto 
the white dwarf is variable during quiescence. 

If the single temperature model is correct, then the im- 
plied accretion rate is in the range 7 x 10~ n to 5 x 10~ 10 
M© yr _1 for distances of 1.6 and 2.45 kpc, respectively. 
This is clearly at odds with the accretion rate required 
to power an outburst in the system every 20 years. If on 
the other hand the cooling flow model is correct, then the 
accretion rate is in the range 2 x 10~ 9 to 1.2 x 10~ 8 M© 
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yr -1 . Allowing for the presence of an absorbed optically 
thick component of the boundary layer that is consis- 
tent with the observations, the accretion rate can be as 
high as 3.6 x 10~ 8 M yr^ 1 . Thus by accounting for 
the more complex nature of the X-ray emission, we have 
shown that the long standing discrepancy between the 
high accretion rate required to power frequent outbursts 
in RS Oph, and the low observed X-ray flux, can in fact 
be resolved. 
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